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Abstract 

The ultra-luminous {Lx ^ lO^^erg/s) intermediate-mass black-hole system HLX-1 in the ESO 243- 
49 galaxy exhibits variability with a possible recurrence time of a few hundred days. Finding the origin 
of this variability would constrain the still largely unknown properties of this extraordinary object. 
Since it exhibits a hardness-intensity behavior characteristic of black-hole X-ray transients, we have 
analyzed the variability of HLX-1 in the framework of the disk instability model that explains outbursts 
of such systems. We find that the long-term variability of HLX-1 is unlikely to be explained by a 
model in which outbursts are triggered by thermal- viscous instabilities in an accretion disc. Possible 
alternatives include the instability in a radiation-pressure dominated disk but we argue that a more 
likely explanation is a modulated mass-transfer due to tidal stripping of a star in an eccentric orbit 
around the intermediate-mass black hole. We consider an evolutionary scenario leading to the creation 
of such a system and estimate the probability of its observation. We conclude, using a simplified 
dynamical model of the post-collapse cluster, that no more than 1/100 to 1/10 of M,< lO** Mq IMBHs 
- formed by run-away stellar mergers in the dense collapsed cores of young clusters - could have a 
fewx 1 Mq Main-Sequence star evolve to an AGB on an orbit eccentric enough for mass transfer at 
periapse, while avoiding collisional destruction or being scattered into the IMBH by 2-body encounters. 
The finite but low probability of this configuration is consistent with the uniqueness of HLX-1. We 
note, however, that the actual response of a standard accretion disk to bursts of mass transfer may 
be too slow to explain the observations unless the orbit is close to parabolic (and hence even rarer) 
and/or additional heating, presumably linked to the highly time-dependent gravitational potential, 
are invoked. 

Subject headings: X-rays: individual(ESO 243-49 HLX-1) - accretion, accretion discs - instabilities - 
stars: binaries: close - galaxies: star clusters - stellar dynamics 



1. INTRODUCTION 

HLX-1 is the b rightest ultra-luminous X-ray (ULX; 
see iRobertsI 120071 for a review) source known, located 
in the outskirts of the edge-on SOa spiral galaxy ESO 
243-49 with a maximum luminosity of '--^10'*^ erg s~^ 
(jFarrell et al.l[2009t iGodet et all 120091 ). The recent dis- 
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covery (jWiersema et al.ll2010f ) in the optical HLX-1 spec- 
trum of an emission line consistent with Ha at the 
redshift of ESO 243-49 (z = 0.0223) irrevocably con- 
firms its association with this galaxy at a distance of 95 
Mpc. With observed X-ray luminosities reaching above 
2o42 ergs~^ HLX-1 is super-Eddington if the black- hole's 
mass is less than ~ lO^M f:^. Beaming effects (e.g. iKingl 
120081 iKording et al.l[200^ have been proposed as viable 
mechanisms for producing the apparent super-Eddington 
luminosities seen from other ULXs. However, beaming 
is unlikely to explain HLX-1 's extreme luminosity due to 
the observed large-scale variability (which appears sim- 
ilar to that seen from Galactic stellar mass black hole 
binaries that are not viewed down the jet-axis) and the 
luminosity of the line (which is an order of magni- 
tude above that expected from reprocessing in the lo- 
cal absorbing material; Wiersema et al. 2010). Based 
on AGN-type scali ng the Hg luminosity might suggest 
a mass < 1500 ModW lersema "eFaniMol ). but since it is 
not clear how the line is related to the accreting system 
this estimate is highly uncertain. By taking the con- 
servative assumption that HLX-1 excee ds the Eddingto n 
limit by no more t han a factor of 10 (|Begelm an ■2002), 
iFarrell et al.l ()2009D placed a lower limit on the black hole 
mass of 500 M© . However, IGodet et al.l (|2010f) obtained 
from a disc-blackbody fit a peak luminosity of 1.3 x 10^^ 
ergs-i (August 29, 2010) and a temperature 2.7 x 10^ K. 
Comparing this value to that obtained from the (non- 
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relativistic) formula for the effective temperature at the 
inner edge of an accretion disk around a black hole: 

rin«6xl0« f^^Y\-V^K, (1) 

where L42 is the luminosity in units of 10"*^ ergs^^, r] = 
O.lryo.i is the accretion efficiency (< 0.4), A/4 the black 
hole mass in units of 10 000 Mg and x = c^R/2GM is the 
radius measured in units of the Schwarzschild radius, one 
concludes that M > IO^'Mq. 

This conclusion is supported by more refined meth- 
ods such as fitting the XMM-Newton European Photon 
Imaging Camera (EPIC) spectra with various models 
(KerrBB, KerrD, slim-disc, BHSPEC) which favor a mass 
M - lO^M© (jGodet et aLllMTillDa^. et alJlMl . The 
same is true of the normalization of the diskbb model. 
One should add that the temperature and luminosity 
observed d uring the previous luminosity peak (on the 
2009-08-16: iGodet et al]|2009D were very close to those 
determined recently. 

After almost two years of monitoring with Swift, one 
sees that the X-ray light-curve of HLX-l shows vari- 
ability with a cha racteristic time of ~ 10^ s (Fig. [TJ 
IGodet et al.l (j2010f )). The available data seem to indicate 
a FRED (Fast Rise Exponential Decay) - type shape. 
The recurrence time seems to be ^ 380 day^ (consis- 
tent also wi th XMMl and two Rosat observations, see 
IWebb et al.l 1201 11 and the decay-time from (the first) 
maximum is ^ 90 days. The rise is very steep and oc- 
curs over a timescale of about one week. The amplitude 
of the X-ray variations is a factor of ^ 20 — 50. The 
decay-from-maximum light-curve shows small "re-flares" 
not dissimilar to those sometimes observed in low-mass 
X-ray transient sources (e.g. Chen et al. 1997). The end 
of the second outburst's light-curve is typical of X-ray 
transients suggesting the propagation of a cooling front. 
In addition, the hardness of the X-ray spectrum of HLX- 
1 follows the trend observed in X-transients: it increases 
with declining luminosity. A recent deep observation 
with XMM-Newton on 2010 May 14 (MJD 55330) con- 
firms that the spectrum does harden at lower luminosities 
(^3 X 10''° erg s~'; Servillat et al. in preparation). The 
analysis of these data suggests that there may be some 
unresolved contribution from the nucleus of ESO 243-49 
in the X-ray emission. However, spectral modeling (us- 
ing a thermal plasma model to represent emission from 
the galaxy) indicates that this is likely to contribute no 
more than ~20% of the total observed flux, and thus that 
the low state luminosity is dominated by emission from 
HLX-l. 

These properties of the HLX-l variability suggest that 
it could result from the same thermal-viscous instability 
that drives the outbursts of X-ray transient systems and 
dwarf-nova stars, and that the HLX-l behaviour might 
be described by a variant of the corre sponding disk in- 
stability model (DIM; see lLasotall200lL for a review). 

2. ACCRETION DISK INSTABILITIES 

2.1. Accretion disk size according to the DIM 

The luminosity of HLX-l (- lO''^ erg s^^) and the 
variability timescale (~ 10^ s) imply a huge accretion 

11 373 days according to lKonel i20Wt 



rate onto the compact object, of the order of lO^'' Mq 
yr^'. Sustaining such a high mass transfer rate excludes 
wind accretion. We will therefore assume that the pu- 
tative stellar companion of HLX-l fills its Roche-lobe, if 
only during part of its orbit, losing matter that forms an 
accretion disk around the black hole. This assumption is 
also supported by the observed intensity related spectral 
changes. The black hole mass is assumed to be around 
10** M0 so the mass ratio q = M^/Mi,, where M, and 
Mi, are respectively the black-hole and stellar-companion 
masses, will typically be q w lO^'' — 10^'^ ^ 1 and the 
Roche-lobe radius can be written 

RL^a{qlif'\ (2) 

where a is the orbital separation (jPaczvnskil 119771 for 
the moment we assume a circular orbit). When q is so 
small, matter circularizes very close to the donor star 
and the 2:1 Lindblad resonance appears at a radius 
0.63a within the primary's Roche lobe (Lin & Papaloizo^ 
Il979t ). a ffecting the accretion d isk formation and struc- 
ture fsee lYungelson et al.l 120061 for a discussion). In the 
following we will assume the disk forms and has a size 
i?D ~ o- 

Here, we assume the variability in HLX-l is related to 
the disk instability model. This requires that the disk is 
big enough to allow its temperature to fall below the tem- 
perature where hydrogen recombines (jLasotall2001|) . Ac- 
cording to the DIM the accretion rate at outburst max- 
imum is nearly constant through the disk and roughly 
equal to the (upper) critical accretion rate (|Lasota et al.l 
1008): 

M+ 1 « 3 X 10-6 a-^-^^Rli^ M^"-'^ M^yv-' , (3) 

where R13 is the disk radius in units of 10"'^'^ cm and 
M, = A/4 X 10"* M0 is the mass of the black hole. This 
limit is essentially independent of the Shakura-Sunyaev 
disk-viscosity parameter a. Therefore the maximum 
(bolometric) luminosity is equal to 

imax « 1.6 X 10^° 770.1 Rft ergs-1 . (4) 

Hence for imax the disk radius will be equal to 

Rf^^"" w 5.2 X 10^^r]a,i^^M^-^'^ cm (5) 

Strictly speaking, for (inside-out) outbursts starting in 
the inner disk regions the radius in Eq. Q corresponds 
to the distance reached by the outside-propagating heat- 
ing front. This distance can be shorter than the actual 
disk outer radius. In such case, the value given in Eq. ([5]) 
is only a lower limit for the disk radius. Note also that 
this assumes a non-irradiated disk. Taking into account 
irradiation of the outer regions of the disk by X-rays from 
the i nner regions always moves the critical radius further 
out (jDubus et al.l Il999f ) . Alternatively, an irradiation- 
dominated disk reprocessing passively C ~ 0.1 of the X- 
ray luminosity would still see its temperature fall below 
8000 K (hydrogen ionization) only for radii greater than 
2x10^'^ cm, similar to the above estimate. So, again, Eq. 
^ is only a lower limit on the disk radius. 
From the Kepler's law a is 

a = 6.3 X 10^^ mI^^P^^^ cm , (6) 
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Figure 1. The 2008 - 2011 Swift-XKT Photon Counting grade 0-12 light-curve of HLX-1 in the 0.3-10 keV energy range. 



with Pd being the orbital period (in days). Combining 
with Eq.Q, the mean density of the companion p can 
be written in terms of the orbital period only 

0.057 P^^gcm-^ (7) 

A disk size i?D ^ 5 x 10^^ cm implies an orbital pe- 
riod > 23 days. The secondary star mean density for 
the hypothetical companion of HLX-1 gives 10"'* g cm"'^ 
suggesting a red giant or a massive supergiant. A pri- 
ori, this is consisten t with severa l possible evolutionary 
scenarios for HLX-1 (jFarrell et a l. 2010, but see Sect. 3). 

2.2. Outburst properties 

The presumed outburst of HLX-1 does not look like a 
"standard" full-blown outburst of a black-hole transient 
low-mass X-ray binary. In all such phenomena, after the 
outburst, the system declines (sometimes after one or 
two re-bounds) to a quiescent, very-low luminosity state. 
According to the DIM the (non-constant) accretion rate 
is then everywhere lower than the critical value given 
by Eq. ([3|) (in fact it is lower than the lower critical 
value M~^^). Assuming that the disk terminates at the 
ISCO (Innermost Stable Circular Orbit) this implies a 
luminosity of < 2 x 10^° erg s~^. Even allowing for inner 
disk evaporation one can only incre ase this luminosity 
by 4 or 5. say, orders of magnitude (iLasota et al.l I1996I -. 
iMenou et all boOOt iDubus et all |200lD, still well below 
the observed ^ 2.6 x 10'*° erg s~*. 

2.3. Irradiation dominated disk 

The moderate amplitude of the outburst in HLX-1 sug- 
gests that only part of the disk is involved i. e. that the 
DIM heating and cooling fronts propagate only in a re- 
stricted domain of the disk. Cooling fronts can be stalled 
when X-ray irradiation from the inner regions of the disk 



prevents cooling in the outer regions. If the disk irradi- 
ation is directly tied to the mass accretion rate onto the 
compact object this does not prevent the disk from emp- 
tying b ut the decay ligh t curve falls dow n exponentially 
(jKing fc Ritter 1998; Dubus et all 1200 More compli- 
cated behavior can arise if X-ray irradiation is not di- 
rectly tied to the inner mass accretion rate. 

Here, irradiation by a hot supergiant star located just 
outside of the disk can have an impact on its stability, in- 
dependently of the mass accretion ra te . Such a possibil- 
ity has been advocated bv lRevnivtsev et al.l ([2002) in the 
case of outbursts of the (presumably) super-Eddington 
source V4641 Sgr. Such irradiation would be constant in 
time (or at least accretion-independent) and its effects 
different from the accretion self-irradiation. Constant ir- 
radiation of an unstable disk typically leads to long term, 
low-amplitude modulation of the inner mass ac cretion 
rate. An example calculation is shown in iDubui ([2005). 
The cooling front propagates inward down to the radius 
where constant irradiation keeps the disk hot. The re- 
gion participating in the instability is limited with little 
of the disk mass accreted at each cycle. 

Although determining the exact timescale and ampli- 
tude of the modulation requires numerical calculations, 
the typical outburst timescale will be linked to the vis- 
cous timescale of the minimum unstable radius (where 
the temperature must be sa 10'* K) 

^ _ ~ 115 a-^T4^*i?Jf A/4' Vears, (8) 

T4 is the disk temperature in units of 10* K and we have 
used the Shakura-Sunyaev prescription for the kinematic 
viscosity coeff icient u = ac ^/il.k- For comparison, the 
case shown in iDubusI (|2005D had a modulation period of 
about 0.5 years but had R13 0.01 and M4 sa 0.001. 
The model appears difficult to reconcile with the fast 
timescale of the variations in HLX-1. The disk must 
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be at least as large as required by Eq. (0) (actually 
larger by a factor 2 when including irradiation heating) 
in order to be unstable in the first place. In order to 
obtain small amplitude outbursts, the radius to which 
the cooling front propagates (Eq. [8]) cannot be orders-of- 
magnitude lower than the outer disk radius. The modu- 
lation timescale will be too long for any reasonable set of 
parameters. We conclude that the variability in HLX-1 
is unlikely to be related to the DIM in the sense that it 
is difficult to reconcile the timescales with a disk large 
enough that its temperature at the outer radius is below 
« 10 K. Therefore the accretion disk in HLX-1 is very 
likely to be hot and thermally stable. 

2.4. Radiation-pressure dominated disk 

We note that in considering the origin of HLX-l's 
variability one should take into account that it is a 
near-Eddington X-ray source. For L/L^dd ^ 0-01 
the standard Shakura-Sunyaev disk is radiation-pressure 
dominated and opacities are due mainly to electron- 
scattering. Such a disk is thermally and viscously 
unstable and the resulting variability has been mod- 
eled producin g FRED-like outbursts with amplitude fac- 
tors < 100 (iTaam fc Lid Il98l iLasota &: Pe"0 ll99lL 
iSzuszkiewicz fc Milleril200l . thus interesting in the con- 
text of HLX-1. The unstable region is limited to the inner 
^ lOO Rg from the black hole s o it is generally short al- 
though [M^eF^^^ringli ()2006f ) find a recurrence time of 
w 250 years for the case of a 10^ black hole accreting 
at 10% of the Eddington rate. Contrary to the case of the 
dwarf-nova DIM, the physics of the radia tion-pressure 
therm al instability is rather controversial (Hiros e et al.l 
[200930). It remains to be seen whether the amplitude 
and recurrence time of the variability in HLX-1 can be re- 
produced by a radiati on-pressure dominated disk around 
a 10^ Mq black hole (|Zheng et al.ll2011l: iCiesielski eFall 
l2OT0h . 

3. ORBITAL ORIGIN OF THE LONG-TERM VARIABILITY 

We have shown above that the disk-instability origin of 
the observed long term variability is very unlikely. Here, 
prompted by the recurrent nature of the X-ray flux varia- 
tions, we consider whether the variability could be orbital 
in origin. 

3.1. Orbital modulation 

Other ULX systems, such as X41.4+60 in M82 and 
NGC 5408 X-1 ar e variable on timescales of 62 and 11 5 
days respectively (jKaaret et al.ll2006l ; lStrohmaveiil2009D . 
These variations have been interpreted as reflecting or- 
bital modulations. As with those systems, the compan- 
ion star of HLX-1 would have to be a massive supergiant 
star (p « 10~^gcm~^) if the 380 day timescale is the 
orbital period. One mechanism that would cause the 
X-ray modulation is absorption and scattering by the 
companion's stellar wind. However, this explanation is 
unlikely for HLX-1 for three reasons. First, the X-ray 
flux is modulated only by a factor ^ 2 in the case of 
the aforementioned ULX sources compared to a factor 
20 - 50 in HLX-1. Second, their light-curves are quasi- 
sinusoidal in contrast with the FRED-like shape of HLX- 
1. Finally, the observations of HLX-1 do not show a sig- 
nificant change in the X-ray absorption column density 



at different luminosity levels (but this is challenging to 
constrain within the statistical quality of the X-ray data 
available) . 

Recently iKingI (|2011h proposed that ULXs in glob- 
ular clusters could be Ultra Compact X-ray Binaries 
(UCXBs) in which a neutron star accretes, at slightly 
super-Eddington rates, matt er los t by a Roche-lobe fill- 
ing companion white dwarf. iKingj (|201lD suggested that 
a version of this model might apply to HLX-1 but left 
the 380-day variability unexplained. In such a scenario 
one could be tempted to identify this variability with the 
so-called super-orbital period observed e.g. in the UCXB 
4U 1820-303, where a variability 170 days is observed in 
addi tion to the orbital perio d of 11 minutes is observed 
(see iZdziarski et al.l l2007bl and references therein) . In 
this system the intrinsic luminosity varies by a factor 
of > 2 only but larger amplitudes could be envisaged 
in the framework of a hierarchical triple system model 
(jZdziarski et al.ll2007a[ ). However, the hardness- intensity 
relation in 4U 1820-303 shows a pattern completely dif- 
ferent from jhat obser ved in HLX-1. Of course, as men- 
tioned by Kin3 pOllf ). if the Ha line is emitted by the 
accretion flow the UCXB model is ruled out. Clearly 
more observations in the optical domain are required to 
decide the viability of this (and other) models. 

3.2. Modulated mass transfer from a donor in an 
eccentric orbit 

A scenario in which the disk in an eccentric binary 
would exist semi-permanently (the viscous time of a 
standard accretion disk being long, see Eq. ([8))) and 
outbursts would be triggered by increased mass-transfer 
rate during the passage of the companion appears to 
be the only serious possibility left. It would be, in a 
sense, the equivalent of the so-called mass-transfer insta- 
bility model advocate d previously for dwarf-nova stars 
(jBath fc Pringld[T98l . In the case of HLX-1 one has to 
show that such an binary system can be formed and sur- 
vive long enough to be observed. Here, we present the 
outline of such a scenario, which we explore in m ore de- 
tail elsewhere (jAlexander. Dubus fc Lasotall201ll ). One 
uncertainty is how eccentric the orbit must be to allow 
fast enough accretion that can reflect the orbital mod- 
ulation. The simple assumptions we make here are not 
expected to hold for very high eccentricity, but they can 
provide an upper limit on the probability of an eccentric 
binary donor in HLX-1. 

3.2.1. The set-up 

For the moment observations provide only very gen- 
eral constraints on possible evolutionary models of our 
hypothetical binary system. As an example we assume 
here that the P^, sa 380 d period observed in the X-ray 
light curve of HLX-1 reflects the orbital period of an 
evolved donor of mass M^, and radius i?^ , in an eccentric 
orbit around an IMBH of M, - 0(10* Mq). To feed the 
accretion (quasi)-periodically, the orbit must graze the 
tidal disruption radius rt ~ i?*(M,/M*)^/'^ at periapse, 
= a*(l — e^), where a* is the orbital semi-major axis 
(sma), and the eccentricity. Therefore one obtains 
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(22.1 AU = 7 X lO^-^cm). The extremely high mass- 
loss requirement of 0(10^* MQyr^^) strongly suggests 
an AGB d onor (initial mass range of ~ 0.5 Mq to 10 M©; 
IHabing fc O lofsson 2003,), which can reach ~ 10 — 
10~^ M (D yr~^ m ass-loss r ate ev en without "tidal induce- 
ment" (jBowen fc Willsonlll991[ ). We will adopt here as a 
fiducial AGB progenitor a M* = 4 Mq star with a main- 
sequence (MS) radius of 2.3 R© and therefore a MS life- 
time of ti, ~ 200 Myr, which expands to R^, ^ 100 R© in 
the AGB phase. The eccentricity required to graze the 
tidal radius is then 



1 - et^rt/ai. 



tO.3 



100 R, 



4Ma 



-1/3 



P. 



380 day 



-2/3 

(10) 



While still on the MS, the donor orbit is stable to both 
tidal and gravitational wave decays, since ~ 75rt ~ 
lO^rg. However, once it reaches the AGB phase, it will 
be quickly destroyed within 10^ — 10^ yr, due to a com- 
bination of the high mass loss rate and the tidally driven 
orbital decay. 

The scenario propose d he r e is d ifferent from the one 
suggested by Hopnian e t al.l ()2004f ) to explain the non- 
periodic ULX in the young cluster MGG-11 in M82. 
There, it was assumed that the progenitor was tidally 
captured while still on the MS, and that the orbit de- 
cayed and circularized by the tidal interaction before the 
star reached the AGB, resulting in a relatively steady 
mass transfer rate on a very short-period circular orbit 
evolving through emission of gravitational radiation. 

3.2.2. Evolution and survival 

Run-away stellar mergers in the dense collapsed cores 
of young clusters are widely considered to b e a natural 



pathway to the formation o f an IMBH (Ebisuzaki et al 
2001; Portegics Zwart & M cMillanI I2002I ; iGtirkan et al 



2004 .Freitag et al., ■2006b.al ). The scenario considered 



here is based on the simulations bv lGiirkan et al.l (j2004 
hereafter GFR04J3- They find that for a wide range 
of initial parameters, a dense cluster with a broad Ini- 
tial Mass Function (IMF) will undergo mass segregation 
and core collapse after a time tec ~ 0-1^^,1/2(0), where 
^r. 1/2(0) is the initial 2-body relaxation at the half-mass 
radius of the cluster. When the core-collapse time is 
shorter than the minimal MS lifetime for massive stars 



tec < mini* ~ 3 - 4 Myr, 



(11) 



the massive stars that segregate in the dense collapsed 
core can undergo run-away mergers on a timescale 
shorter than the stellar evolutionary timescale. If the 
merged "super-star" can cool fast enough and avoid rapid 
mass loss, it may collapse to form an IMBH. 

Guided by the GFR04 simulations, in what follows we 
will assume that this process forms an IMBH with a mass 
ratio Mt/Mc — 0.002, where Mc is the total initial mass 

A scenario in which a nucleated dwarf galaxy is undergoing a 
period of accretion due to a recent passage through the host galaxy 
is often mentioned in the context of HLX-1. Although the presence 
of IMB H in nuclei of dwarf galaxies has still to be demonstrated 
(see e.e. lReines et aLl l2011) we think we are justified in considering 
this possibility just as a slightly different realization of the scenario 
we are studying. 



of the cluster. As reference we will use Model 2 of GFR04 
in which a Salpeter IMF extending from 0.2-120 M© is 
assumed (mean mass {Mi,) = 0.7 M©, maxM^,/ {M^) = 
174). The model assumes a Plummer distribution with 
a length-scale = (2^/3 - 1)1/2^/3 ~ 0.766ri/2, where 
^1/2(0) is the initial half mass radius. The core-collapse 
time is: 

Nk 1 

log lO-^iVe Vpi 



where po = (10^ M© pc '^)pg is the initial central density 
and Nc = lO^A^g is the stellar number. The condition 
tec < 3 Myr then translates to po > 10^ M© pc~'^ (for 
Mc = 5 X lO^M©, Nc = 7.25 x 10^, R^ = 0.11 pc). Al- 
though young clusters and super star clusters with relax- 
ation times well belo w 30 Myr are observed (jFiger et al.l 
l2002t[H"o fc Filippen ko 1996) it is less clear whether such 
extreme central densities are realized in nature. A lower 
IMBH mass will ease this problem since according to Eq. 
(HH), Po ~ Ml 

In such dense systems mass segregation proceeds very 
rapidly and accelerates the evolution toward core col- 
lapse. Following the results of GFR04 (their Fig. 5) 
we take for the stars closest to the newly formed IMBH 
(M*)^ = 4M©, with a corresponding radius (i?*)^ = 
2.3411©, as the mean stellar mass and radius within 
the IMBH radius of influence, — GM^/a^, where 
is the typical ID velocity dispersion of the host clus- 
ter (cr2(o) = GMc/QRc ^ 185kms-i for the Plummer 
model here). It is noteworthy that mass segregation ini- 
tially concentrates intermediate mass AGB progenitors 
near the IMBH. 

Following the collapse the dynamical response of the 
cluster by mass segregation and expansion leads to 
the formation of a relaxed stellar cusp inside r/i, with 
nnir) oc r~"" {aa = 7/4-5/2) for the massive stars, and 
nijr) oc r~^^^ for the low-mass s t ars (iBahcall fc Wolj 
1977t [A lexander fc HopmanI l2009t iKeshet et all 120091 ; 

101). 



Preto fc Amaro-Seoane 20101). The cusp extends inward 
down to the collision radius, ri„ = max(rt, rcoii) 



'coll 



(Af,/Af*)i?*, at which 2-body relaxation ceases to be ef- 
fective because velocities are so high that only physi- 
cal collisions can sub stantially change the stellar orbits 
(|Frank fc Reed fT976[ ). For a standard MS mass-radius 
relation this gives 



rcoii ~ 48.7(Af./104M©)(M,/Af©) 



-0.43 



AU, 



(13) 



so that for the cusp to extend down to a* ~ 22 AU the 
stars must be more massive than ~ 3.7 Af©. One notices 
that this lower limit coincides with the mass range pre- 
dicted by pre-IMBH mass segregation, while being low 
enough to include AGB progenitors. 

Equilibrium is established when the flux of gravita- 
tional binding energy that is released when stars are 
destroyed by the IMB H equals the flux c arried by the 
expanded cluster core (|Heggie et al.l [20071) . Assuming a 
Plummer initial distribution one obtains the mass of the 
cusp Mh as (e.g. iBaumgardt et al.ll2004l ) 



Mh ~ [648/(3 - aH)]iM,/McfM, 



(14) 



In contrast to the radius of influence of a supermas- 
sive black hole (SMBH) in an approximately isothermal 
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galactic nucleus, that of an IMBH contains only a very 
small number of stars, Nu = Mh/ (M*)^ ^ 0(10). 

Equations ([12]) and ((T4)) , together with the system pa- 
rameters Af., M,/Mc, aH, {AQ and J fully 
describe the IMBH cusp structure and properties in this 
simplified model, and allow calculating the mean num- 
ber of stars on orbits with a < and e > e^,, {Na,e) = 
Nh[{5/4:)ai,/rfi]'^~°'" {1 — ej), where an isotropic cusp is 
assumed. The Poisson probability for having at least one 
star on a donor orbit is then Pi = 1 — exp[— (iV(j,e)]- 

The cusp mass and its stellar density rise with the 
IMBH mass (Eq. HM . As a consequence. Pi increases 
with A/,, but with it also the rate of destructive stel- 
lar collisions, a nd the drain rate (star- star scatterings 
into the IMBH, [Alexander fc LiviolboO^ . O ver the rel- 
evant range of IMBH masses, Pi rises from < 10"'^ to 
> 0.1. The IMBH mass that maximizes Pi subject to 
both the collisional and drain constraints lies in the range 
of few X IO^Mq, with max Pi ^ few x 0.01. 

The observed residual UV fluxes from HLX-1, af- 
ter subtracting a disk model and correcting for ex- 
tinction (Eb-v = 0.042), are 4.8 x 10"^^ ergcm"^ s^^ 
{Lnijv — 5.2 x lO^^ergs"^ for isotropic emission at 
D = 95Mpc) in the NUV (2147 A - 3467 A), and 
4.7 X 10-15 erg cm-2s-i (Lpuv = 5.1 x lO^^eigs"^ ) in 
the FUV (1233 A - 1821 A) (based on a preliminary anal- 
ysis of HST data, Farrell et al., in preparation). These 
place constraints on the properties of the hypothesized 
birth cluster of the IMBH. 

Stellar population synthesis models together with 
model atmospheres allow to predict the UV flux as func- 
tion of the IMF, metallicity and age of the system. Gen- 
erally, the older the cluster, the less UV it emits. An 
older cluster directly implies a lower mass AGB pro- 
genitor, that can spend longer on the MS. Preliminary 
modeling, both assuming black body spectra, and us- 
ing detailed ste llar atmosphere models (the STARS code, 
lSternber3ll998i) indicates that a minimal cluster age of 
~ (0.3 — 0.6) Gyr yr is required for the cluster's UV lu- 
minosity to fall below the residuals (Assuming a cluster 
mass of 5 X 10^ M© and solar metallicity). This cor- 
responds to an AGB progenitor of ~ (2.7 — 3.5) M©. 
Such longer-lived progenitors are more susceptible to col- 
lisional destruction and scattering into the IMBH, but 
are still probable at the Pi ~ few x 0.01 level. The con- 
straints on the AGB progenitor mass could be relaxed 
somewhat by assuming a higher value of AI,/Mc, and 
possibly by assuming a different metallicity. This re- 
quires a more systematic study. 

After mass transfer starts in the AGB phase, 2-body 
perturbations from other stars can be neglected, since 
the relaxation time at is th ~ few x 10 yr, and the 
timescale to significantly affect the orbital eccentricity, 
(1 — el)th, is longer than the ^ 10"* yr maximal lifetime 
of the donor, for all relevant values of the eccentricity. 

Thus, while an IMBH with Ah 10** M© is not ex- 
cluded in the context of the eccentric donor scenario, 
observations and theoretical considerations favor a some- 
what less massive IMBH with M, ^ few x 10^ A/,, still 
consistent with the constraints derived from the observed 
accretion emission (see Sec. [T]). 

To summarize, we used a simplified model, based on 
results from iV-body simulations of runaway merger for- 
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Figure 2. A comparison of the observed de-reddened limits of 
the residual UV fluxes (after subtraction of a disk model) to stellar 
population synthesis model predictions for different metallicities 
(Z) and low mass cutoffs (Mi) and assuming a 5 X 10^ M0 cluster, 
as function of the maximal progenitor mass (the turn-off mass), 
which sets the age of the cluster (also shown) and hence its minimal 
UV luminosity. 

mation of a seed IMBH, to predict the stellar distribution 
around a newly formed IMBH. We estimated the prob- 
ability of finding an AGB progenitor on an orbit that 
could explain the long period variability of HLX-1 in 
terms of mass transfer from an eccentric evolved donor. 
We find that between 1/100 to 1/10 of IMBHs with 
10'' M© could have a fewx 1 M© MS star evolve to 
an AGB while avoiding collisional destruction or being 
scattered into the IMBH by 2-body encounters. An im- 
portant caveat is that the validity of the model and our 
conclusions are limited by the neglect of the longer-term 
evolution of the cusp and host cluster, and in partic- 
ular by the neglect of the eff ects of 2-body r elaxation 
and fast resonant relaxatio n (iRauch fc Trem aine 1995 
IHopman fc Alexanderi [2006[ ) on the distribution of the 
stars around the IMBH. 

The low probabilities are fully consistent with the ap- 
parent uniqueness of HLX-1. Theoretically we can only 
make statements about the conditional probability since 
we do not know the space density IMBH clusters but 
observat iona l ly, de spite large-scale searche s (in Chandra 
data by [03 (|2011[ ) and in XMM data by IWalton et all 
(|2010[ )). there have not been any other HLX-1 like ob- 
jects found. Most probably HLX-1 is alone in the local 
Universe. 

3.2.3. Variability timescales 

In our scenario involving a 1 — e = 0.3 binary in which 
a star circles a ^ lO^'M© black hole on an 380 day orbit, 
the periapse will be at ^ 10^^ cm. The mass loss and 
disk formation processes is usually described by the star 
filling an "instantaneous" Roche lobe at p eriapse (see e.g. 
iSepinskv et"aIll2007HLaioie fc Sillsll201lL in the context 
of high mass X-ray binaries). This description is increas- 
ingly inaccurate as the eccentricity grows but matter lost 
during the periapse passage is expected to circularize at 
about this distance. The transferred material will have 
to diffuse inwards, on a viscous timescale. From the dis- 
cussion in Sec[21 the disk is most likely hot with a tem- 
perature > 10, 000 K at the outer edge. The diffusion 
timescale for the transferred material would be hundreds 
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of years rather than days according to Eq. [S] As a re- 
sult, the accretion timescale will smooth out the bursts 
in mass transfer. 

One way out of this difficulty is to increase the disk 
temperature. Irradiation is unlikely to increase the disk 
temperature much: the irradiation temperature Tin- — 
(CLxI'^-Kcr'i^R^Yl^ is w 3500 K at lO^'' cm (where 
C K, 10^'^ parametriz es our igno rance of the irradiation 
geometry and albedo. Ibubus et~ al. 1999). Incoming ma- 
terial may shock heat the outer disk to high tempera- 
tures, especially in the case of a tidal disruption on a 
parabolic orbit when the material ejected from the star 
has a large range of velocities relative to the Keplerian 
disk veloci ty, up to the escape speed of the star (see e.g. 
lRee^[T988h . The response of the disk to a burst of mass 
transfer in this situation has not been modeled. Note, 
however, that in close low-mass binary systems the "hot 
spot" resulting from the mass-transfer stream impact sig- 
nificantly heats the outer disk without drastically affect- 
ing the dynamical properties of the disk because the ther- 
mal timescale on which the extra heating is radiated away 
in a thin disk is much sho r ter than the viscous timescale 
(|Buat-Menard et al.l[200l[Smaa[200l . 

Another way to reduce the diffusion timescale is to 
reduce the periapse distance. The viscous timescale is 
already down to 3 years if periapse is at rt ~ 10^^ cm 
(assuming the temperature varies as T - in Eq. 

[S]). In such a model the disk does not disappear com- 
pletely in between mass transfer episodes (in accord with 
the non-zero minimum flux) and the impulsive increase in 
mass transfer leads to increase in lightcurve on timescales 
more like a fraction of tms (^ 0.1) and then decays on 
tms or so. The passage of the star at periapse may also 
lead to the excitation of waves in t he disk that will en- 
hance angular momentum transport (|Spruitlll987f ) . Tidal 
waves may provide an additional source of heating in the 
disk. The price to pay is that such a binary has a much 
lower probability because (1) the orbit must be nearly 
parabolic, e — ^ 1, (2) the depletion of the phase space 
density of orbits near the loss cone, which was not taken 
into account in our simple estimate (sec. I3.2.2|) . and (3) 
the unstable nature of such an orbit. In fact, this in- 
stability could lead to an observed complete disruption 
within a few orbits (few years). 

4. DISCUSSION AND CONCLUSIONS 

The extremely high luminosity, light curve shape and 
X-ray spectrum evolution of HLX-l point toward disk 
accretion around a 10** Mq black hole fueled by a Roche- 
lobe filling star. We have examined the conditions under 
which the X-ray variability might be explained by the 
disk instability model. We find this requires an accre- 
tion disk much too large for corresponding timescales to 
be compatible with the observed X-ray variability. Any 
accretion disk around HLX-1 is most likely small enough 
to be hot and stable against the DIM. One cannot ex- 
clude that the variability is due to the instability that 
can arise in radiation-pressure dominated disks. How- 
ever, the physics behind this instability is not well known 
and it is not clear whether this will lead to the correct 
outburst amplitude and timescales. 

The variability would be much easier to explain with 
a stellar-mass black hole but the luminosity would ob- 
viously be problematic. Conversely, the luminosity is 



no issue for a supermassive black hole but the ampli- 
tude aiid__Winescales would be an insurmountable prob- 
lem (|Hameurv et al.ll2009[ ). 

We have shown, however, that a viable description of 
the HLX-1 variability can be provided by a model in 
which enhanced mass transfer into a quasi-permanent 
accretion disk is triggered by the passage at periapse of 
an evolved (AGB) star circling the IMBH on an eccen- 
tric orbit. Using a simplified model based on the results 
of N-body simulations we concluded that such systems, 
although not common, are realistically observable. How- 
ever, the actual response of a standard accretion disk to 
bursts of mass transfer may be too slow to explain the 
observations unless the orbit is close to parabolic and/or 
additional heating, presumably linked to the highly time- 
dependent gravitational potential, is invoked. 

The validity of our conclusions is limited by the ne- 
glect of the longer-term evolution of the cusp and host 
cluster, in particular by neglecting the effects of 2-body 
relaxation and fast resonant relaxation on the distribu- 
tion of the stars around the IMBH. In general, the very 
small number of stars in the IMBH cusp casts doubts on 
the applicability of the statistical approaches commonly 
used to analyze dynamics around SMBHs. In addition, 
the comparatively high density of unbound (cluster) stars 
in the cusp complicates the analogy with known results 
from SMBH cusp dynamics. Progress in the analysis of 
the post-formation evolution of IMBHs will require time- 
dependent modeling and N-body simulations. 

The further evolution of the intriguing variability pat- 
tern of HLX-1 (periodicity, amplitude) as well as obser- 
vations of its optical counterpart should shed light on the 
origin and nature of this extraordinary system. 
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